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Abstract - We show that when a non-wetting fluid drains a stratified porous medium at suffi- 
ciently small capillary numbers Ca, it flows only through the coarsest stratum of the medium; 
by contrast, above a threshold Ca, the non- wetting fluid is also forced laterally, into part of the 
adjacent, finer strata. The spatial extent of this partial invasion increases with Ca. We quantita- 
tively understand this behavior by balancing the stratum-scale viscous pressure driving the flow 
with the capillary pressure required to invade individual pores. Because geological formations 
are frequently stratified, we anticipate that our results will be relevant to a number of important 
applications, including understanding oil migration, preventing groundwater contamination, and 
sub-surface CO2 storage. 
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x/i Introduction. — Drainage, the displacement of a wet- 
# ^ting fluid from a porous medium by an immiscible non- 
wetting fluid, arises in many technological problems, in- 
^cluding oil migration and recovery [lp] , waste CO2 seques- 
^^tration [3j|6], and groundwater contamination [7jj9]. The 
1— ability to accurately predict the flow behavior of the non- 
wetting fluid is critically important in all of these exam- 
^ples. To displace the wetting fluid from a pore, a threshold 
£f~j capillary pressure must build up in the non- wetting fluid 
l/^at the pore entrance; this pressure is given by P 1 ~ 7/at, 
where 7 is the interfacial tension between the fluids and a t 
^^is the radius of the pore entrance 10 -12 . For a homoge- 



rH neous porous medium, characterized by pores of a single 
average size, P 1 is typically much larger than the viscous 
pressure associated with flow into a pore. Consequently, 
. . the flow path taken during drainage depends primarily on 
K^the slight pore-scale variations of a t 13 16 19 . However, 



k^many porous media are stratified, consisting of parallel 



20 22 



$Jj strata characterized by different average pore sizes 

Such additional variation in the pore structure, on scales 
much larger than a single pore, may strongly modify the 
flow behavior 16 18 23 29 . Despite its enormous prac- 



tical importance, a clear picture of how the subtle inter- 
play between capillary and viscous forces determines the 
flow through a three-dimensional (3D) stratified porous 
medium remains elusive. This requires direct visualization 
of the multiphase flow, both at the scale of the individual 
pores and the overall strata. Unfortunately, the medium 
opacity typically precludes such visualization. As a re- 



sult, knowledge of how exactly drainage proceeds within 
a stratified porous medium is missing. 

Here, we use confocal microscopy to investigate drainage 
within a 3D porous medium having parallel strata ori- 
ented along the flow direction; this allows us to directly 
visualize the multiphase flow at pore-scale resolution. We 
find that for sufficiently small capillary numbers, Ca, the 
non-wetting fluid flows only through the coarsest stratum 
of the medium. By contrast, above a threshold Ca, the 
non-wetting fluid is also forced laterally into part of the 
adjacent, finer strata. By balancing the viscous pressure 
driving the flow with the capillary pressure required to in- 
vade a pore, we show how both the threshold Ca and the 
spatial extent of the invasion depend on the pore sizes, 
cross-sectional areas, lengths, and relative positions of the 
strata. Our results thus help elucidate how the path taken 
by the non- wetting fluid is altered by stratification in a 3D 
porous medium. 

Experimental Methodology. — We prepare rigid 
3D porous media by lightly sintering densely-packed hy- 
drophilic glass beads, with polydispersity « 4%, in thin- 
walled rectangular quartz capillaries; these have cross- 
sectional areas A = 1,3, or 4 mm 2 . The packings have 
porosity <\> ~ 0.41, as measured with confocal microscopy 
[l0] . To create stratified porous media, we arrange the 
beads into parallel strata characterized by different bead 
sizes and the same length L; the interface between the 
strata runs along the direction of fluid flow. To enable 



p-1 



S. S. Datta and D. A. Weitz 



Dyed Wetting Fluid Bead 



Area A 




Microscope 
Objective Lens 



warn 



Fig. 1: (Left) Schematic of a stratified porous medium, with 
parallel strata comprised of differently-sized beads. (Right) 
Optical section acquired within a porous medium with a coarse 
and a fine stratum; the pore space is saturated with the 
fluorescent ly-dyed wetting fluid, and the black circles show 
cross sections of the beads making up the medium. 



imaging of the multiphase flow within the media, we use 
fluids whose compositions are carefully chosen to match 
their refractive indices with that of the glass beads; the 
wetting fluid is a mixture of 91.4 wt% dimethyl sulfoxide 
and 8.6 wt% water, while the non-wetting oil is a mixture 
of aromatic and aliphatic hydrocarbons. The viscosities 
of the wetting fluid and the oil are \i w = 2.7 mPa.s and 
\i Q = 16.8 mPa.s, respectively. The interfacial tension be- 
tween the two fluids is 7 = 13.0 mN/m, as measured using 
a du Noiiy ring. We use confocal microscopy to estimate 
the three-phase contact angle between the wetting fluid 
and glass in the presence of the oil, 6 « 5° [To] . 

Prior to each experiment, the porous medium is evacu- 
ated under vacuum and saturated with CO2 gas; this gas is 
soluble in the wetting fluid and prevents the formation of 
any trapped bubbles. We then saturate the medium with 
the wetting fluid; to visualize the fluid, we dye it with fluo- 
rescein. We visualize the pore structure in 3D using a con- 
focal microscope, acquiring a 3D stack of 39 11 /im-thick 
optical slices, each spanning a lateral area of 912 jam x 
912 /im, spaced by 6 jam along the z-direction, within the 
porous medium, away from its boundaries [schematized 
in Fig. 1]. We identify the glass beads by their contrast 
with the dyed wetting fluid; an example slice is shown in 
Fig. 1 [To) . Moreover, we acquire stacks at multiple loca- 
tions along the length of the medium; by combining these 
stacks, we obtain a map of the pore structure of the entire 
medium. 

To investigate the drainage process, we subsequently 
flow the oil at a constant volumetric rate Q, re-acquiring 
3D image stacks at the same positions within, and along 
the length of, the medium. The oil is undyed; we thus 
identify it by its additional contrast with the dyed wetting 
fluid in the measured pore volume. This enables us to 
directly visualize the multiphase flow, both at the scale of 
the individual pores and the overall strata. 
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Fig. 2: Optical sections through an entire porous medium with 
a coarse and a fine stratum, obtained after flow has reached an 
unchanging steady state; imposed flow direction is from left to 
right, (a) Oil flows exclusively through the coarse stratum for 
Ca = 7.2 x 10 -5 ; (b-c) oil also invades part of the fine stratum, 
up to a length Lf from the inlet, for Ca = 1.1 x 10 -4 and 
4.1 x 10 -4 , respectively. The pore space is initially saturated 
with the fluorescent ly-dyed wetting fluid; black circles are the 
beads comprising the medium, while additional black fluid is 
the invading oil. Inlet is at position x = 0.(d) Fluorescence 
intensity in the fine stratum, integrated along y and z, as it 
varies with x for Ca = 1.1 X 10 -4 and 4.1 x 10 -4 as shown 
by the top and bottom traces, respectively. Smooth curves 
show filtered data, and arrows show distance invaded by the 
oil Lf, determined from the inflection points of the curves. The 
two traces, and corresponding curves, are vertically shifted for 
clarity. 



Results and Discussion. — The effect of stratifica- 
tion on the flow behavior is exemplified by draining a 
porous medium with a coarse and a fine stratum, com- 
prised of beads with radius a c = 38 jam and df = 19 /im, 
respectively [Fig. 1]. The medium has length L = 4.2 mm 
and a total cross-sectional area A = 1 mm 2 ; the coarse 
stratum occupies an area A c = 0.44A We flow the oil at 
a constant rate Q = 0.2 mL/hr, corresponding to a capil- 
lary number Ca = ja Q (Q/A)/j = 7.2 x 1(T 5 [30j. The oil 
invades the medium through a series of abrupt bursts into 
the pores, indicating that a threshold capillary pressure 
must build up in the oil before it can invade a pore [To] . 
Interestingly, the oil flows exclusively through the coarse 
stratum, as shown by the optical slice in Fig. 2(a), over 
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Fig. 3: For sufficiently small Ca, oil does not invade into fine 
stratum (open symbols), while for sufficiently large Ca, oil par- 
tially invades the fine stratum as well (filled symbols). Data are 
shown for different porous media, with different cross-sectional 
areas, lengths, and bead sizes; these are thus characterized by 
different values of Ca*, defined in Eq. (1). The threshold for 
invasion into the fine stratum is approximately Ca — 1.7Ca* 
(dashed line). 



an observation time of 30 min [3l] . 

To further explore the drainage process, we increase the 
oil flow rate to Q = 0.3mL/hr, corresponding to Ca = 
1.1 x 10 -4 . Surprisingly, in contrast to the Ca = 7.2 x 
10 -5 case, the oil invades the fine stratum, although only 
partially, over a limited distance, Lf, from the inlet, as 
shown in Fig. 2(b). After an unchanging steady state is 
reached, we flow the oil at even higher flow rates and probe 
the resulting steady state invasion patterns. Interestingly, 
the fine stratum remains only partially invaded for the 
entire range of Ca explored; however, we find that Lf 
increases with increasing Ca, as shown in Fig. 2(b-c). 
This observation contradicts the idea that the fine stratum 



is completely impervious to the oil 32 . 

To quantify the partial oil invasion into the fine stra- 
tum, we integrate the fluorescence intensity in the fine 
stratum along both the y and z-directions, for each posi- 
tion x. Two examples, corresponding to the invasion pat- 
terns shown in Fig. 2(b) and (c), are shown by the upper 
and lower traces in Fig. 2(d), respectively. To determine 
the distance invaded by the oil into the fine stratum, Lf, 
we apply a low-pass filter to these data 38 and determine 



the distance from the inlet to the inflection point of each 
filtered curve [points in Fig. 2(d)]. Consistent with the 
optical slices shown in Fig. 2(b-c), we find that Lf in- 
creases steadily with increasing Ca, as shown in the inset 
to Fig. 5. 

To understand this complex flow behavior, we analyze 
the distribution of pressures in the oil as it displaces the 
wetting fluid. For the oil to invade a pore formed by beads 
of radius a, the capillary pressure at the oil-wetting fluid 
interface must exceed P 1 = 2jcos6/a t , where a t « 0.19a 



Fig. 4: Time sequence of confocal micrographs, taken at 
Ca > Ca*, showing flow of oil (black) from the coarse stra- 
tum laterally into the fine stratum, indicated by arrows in first 
panel; Ca = 4.7 x 10~ 4 and Ca* = 3.8 x 10~ 5 . Imposed flow 
direction is from left to right. 



a pore of the coarse stratum, P 1 ^ c ~ l/a c , is thus smaller 
than that required to invade a pore of the fine stratum, 
P 1 j ~ 1/a/; consequently, we expect the coarse stratum 
to be drained first. This expectation is in direct agreement 
with the observed drainage behavior, shown in Fig. 2(a). 

An additional viscous pressure, P v , drives the continued 
flow of oil through the coarse stratum; we use Darcy's law 
to estimate this as P v (x) — ja (Q/A c )(L — x)/k c [20]. We 
estimate the permeability of the coarse stratum to the oil, 
k c , using the Kozeny-Carman relation k c 



45 Cl- 



10 35 , 37 ; the relative permeability k quantifies the per- 
meability reduction resulting from trapping of the wetting 
fluid within the crevices of the medium, as visible in Fig. 
2(a). We independently measure k « 0.16 using a ho- 
mogeneous porous medium constructed and drained in a 



manner similar to the experiments reported here 10 36 . 

We hypothesize that the oil begins to invade the fine 
stratum when the flow rate is sufficiently large for the vis- 
cous pressure at the inlet, P v (0), to balance the capillary 
pressure required to invade a pore of the fine stratum P 7j /. 
In non-dimensional form, this criterion is 



Ca « Ca* 



2a r cos a r 



0.19a/ L 45 (1 



A 



(1) 



10 12 , 33 34 . The capillary pressure required to invade 



We therefore expect drainage through only the coarse stra- 
tum for Ca < Ca*; for Ca above this threshold, the oil can 
also begin to invade the adjacent fine stratum, in agree- 
ment with our observations [Fig. 2(a-c)]. To test this pre- 
diction quantitatively, we repeat the experiments on many 
different stratified porous media, varying the bead sizes a c 
and ay, medium length L, and the cross-sectional areas A 
and A c \ this enables us to vary Ca* over one order of mag- 
nitude, Ca* ~ 10" 5 - 10" 4 . For all of the media tested, 
we observe exclusive drainage through the coarse stratum 
below a threshold value of Ca, while above this threshold, 
the oil also begins to invade the adjacent fine stratum, as 
shown by the open and filled symbols in Fig. 3, respec- 
tively. We find that the threshold for invasion into the 
fine stratum is given by Ca « 1.7 C a* over a broad range 
of Ca* [Fig. 3, dashed line], in close agreement with our 
prediction [Eq. 1]. 
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Fig. 5: (Inset) Distance invaded into the fine stratum, L/, in- 
creases with increasing Ca; different symbols represent differ- 
ent media characterized by different values of Ca* . Symbols are 
the same as in Figure 3. L is the length of the medium. (Main 
panel) Data collapse when Ca is rescaled by Ca*, and agree 
well with the theoretical prediction Lf/L = 1 — 1.7Ca* /Ca 
(dashed line). 



Within this picture, for sufficiently large Ca, oil is forced 
into the fine stratum not only from the inlet, but also 
laterally, from the adjacent coarse stratum [42H44]. By 
directly visualizing the drainage dynamics at Ca > Ca*, 
we confirm this lateral flow, as indicated by the arrows 
in Fig. 4. We therefore expect that the oil invades the 
fine stratum for all x < Lf, where P v exceeds P 1 j 
Balancing these pressures yields 
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L f /L= 1-Ca*/Ca 



(2) 



We test this prediction by measuring the variation of Lf 
with Ca for the different stratified porous media, char- 
acterized by Ca* ~ 10~ 5 - 10" 4 . For all of the exper- 
iments, we find that Lf/L increases with increasing Ca 
[inset to Fig. 5], consistent with our expectation. More- 
over, the data for different porous media collapse when 
Ca is rescaled by Ca*, as shown in Fig. 5, in agreement 
with Eq. 2. The findings presented in Fig. 3 suggest that 
Ca* should be replaced by 1.7 C a* in Eq. 2; this yields an 
excellent fit to the data, as shown by the dashed line in 
Fig. 5. Our model thus captures both the onset and the 
spatial extent of oil invasion into the fine stratum. 

To test the generality of our results, we also study media 
with three different strata: a coarse stratum, comprised of 
beads with radius a c = 75/im, a fine stratum, comprised of 
beads with radius a/ = 19/im, and an intermediate stra- 
tum separating the two, comprised of beads with radius 
a m = 38/im. We observe flow behavior similar to the case 
of two strata: for Ca = 1.8 x 10" 5 and 4.8 x 10" 5 , the oil 
flows through the entire coarse stratum, and also partially 
invades the intermediate stratum [45]; the spatial extent 
of this invasion increases with increasing Ca. Moreover, 
at an even higher Ca = 9.6 x 10 -5 , the oil also partially 




Fig. 6: Optical section through part of a porous medium with 
three strata, obtained after the flow has reached an unchang- 
ing steady state. Both coarse and intermediate strata are com- 
pletely invaded by the oil (black) , while the fine stratum is only 
partially invaded, as indicated by the arrow. Imposed flow di- 
rection is from left to right; Ca — 9.0 x 10 -5 . The threshold 
Ca* = 2.4 x 10~ 5 and 5.5 x 10~ 5 for the intermediate and fine 
strata, respectively. 



invades the fine stratum. The partial invasion into the in- 
termediate stratum requires the lateral flow of oil from the 
coarse stratum; we thus expect that, if the positions of the 
intermediate and fine strata are switched, the intermediate 
stratum becomes completely, not partially, invaded above 
a threshold Ca. To test this idea, we study a medium 
with the fine stratum separating the coarse and the inter- 
mediate strata. At a Ca = 4.5 x 10 -5 , the oil completely 
invades both coarse and intermediate strata [45], in con- 
trast with the previous case, and in direct agreement with 
our expectation [top and bottom strata in Fig. 6]. More- 
over, at an even higher Ca = 9.0 x 10 -5 , the oil partially 
invades the fine stratum [arrow in Fig. 6] , consistent with 
the picture presented here. These observations confirm 
that the flow path taken by the oil depends not only on 
the geometry of the individual strata, but also on their 
relative positions. 

Conclusions. — Using direct visualization by confo- 
cal microscopy, we demonstrate how stratification alters 
the path taken by a non-wetting fluid as it drains a 3D 
porous medium. For sufficiently small Ca, drainage pro- 
ceeds only through the coarsest stratum of the medium; 
above a threshold Ca, the non- wetting fluid is also forced 
laterally, into part of the adjacent, finer strata. Our results 
highlight the essential role played by pore-scale capillary 
forces, which are frequently neglected from stratum-scale 
models of flow, in determining this behavior. Because geo- 
logical formations are frequently stratified, we expect that 
our work will be relevant to a number of important appli- 
cations, including understanding oil migration 46|47 , pre- 
venting groundwater contamination f7p] , and sub-surface 
storage of C0 2 [48]. 
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